APS/123-QED 



in 
o 
o 

(N 

a 



(N ■ 

> : 
i> 
\o ■ 
o ■ 

O : 
in 

o ■ 

:^ 

Oh! 



A comparative analysis of tunneling time concepts: Where do transmitted particles 

start from, on the average? 

N. L. ChuprikoAQ 

Tomsk State Pedagogical University, 634041, Tomsk, Russi^ 
(Dated: February 1, 2008) 

In this paper we compare the concept of the tunneUng time introduced in quant-ph/0405028 with 
those of the phase and dwell times. As is shown, unlike the latter our definition of the transmission 
time coincides, in the limit of weak scattering potentials, with that for a free particle. This is valid 
for all values of the particle's momentum, including the case of however slow particles. All three 
times are also considered for a resonant tunneling. In all the cases the main feature to distinguish 
our concept from others is that the average starting point of transmitted (reflected) particles does 
not coincide with that of all particles. One has to stress here that there is no such an experiment 
which would give coordinates of all the three points, simultaneously. For measuring the position of 
the average starting point of transmitted particles we propose an experimental scheme based on the 
Larmor precession effect. 

PACS numbers: 03.65.Ca, 03.65.Xp 



I. INTRODUCTION 

As is known [1 II H S IE H , the question of the time 
spent by a quantum particle in the barrier region is still 
controversial. At present there are many different defini- 
tions of this quantity, however none of them is commonly 
accepted. At the same time, for a given potential and ini- 
tial state of a particle, i.e., in the standard setting of the 
tunneling problem, the above question is evident to imply 
an unique answer. 

In our previous paper we have introduced the con- 
cepts of the transmission and reflection times based on 
a separate description of the transmission and reflection 
processes. By our approach, tunneling is a combined ran- 
dom process to consist from two alternative elementary 
ones, transmission and reflection. For a given potential 
and initial state of a particle, we have found an unique 
pair of solutions to the Schrodinger equation, which de- 
scribe separately these processes. This permits one to 
follow the centers of "mass" (CMs) of the transmitted 
and reflected wave packets, at all stages of scattering, 
and, as a consequence, to calculate the times spent by the 
CMs in the barrier region. These characteristic times are 
treated in our formalism as the (average) transmission 
and reflection times for a particle. 

Note, unlike the standard wave-packet analysis 
(SWPA) our approach predicts that the average starting 
point of transmitted (and reflected) particles does not co- 
incide with the starting point of the CM of the incident 
wave packet to describe all particles. By this reason, our 
transmission and reflection times differ essentially from 
the corresponding phase times derived in the SWPA. 

We have to stress that our approach agrees entirely 
with the foundations of quantum theory, for the wave 
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functions for transmission and reflection to underlay it 
are solutions to the Schrodinger equation. Besides, one 
has to bear in mind that the average starting point of 
transmitted particles and that of all scattering particles 
cannot be measured simultaneously, in principle. 

In this paper we show that the formal comparison of 
the well-known phase and dwell times with ours speaks 
in favor of the latter. And, what is more important, 
our formalism can be experimentally verifled. Due to 
the Larmor precession of the particle's spin in an exter- 
nal magnetic field applied beyond the barrier region, one 
can measure, in principle, the average starting point of 
transmitted (or reflected) particles. 

Note, apart from the phase time, the dwell and Lar- 
mor time are, perhaps, the most cited concepts of the 
tunneling time. Thus, with taking into account that the 
Larmor time coincides with the dwell time, it is useful to 
compare our concept of the tunneling time with those of 
the phase and dwell times. 



II. CRITIQUE OF THE DWELL-TIME 
CONCEPT 

In we pointed to the principal shortcomings of the 
phase time concept. Now, before a formal comparing of 
the three concepts, we want to dwell on the principal 
shortcomings of the dwell time concept. 

We begin with the fact that calculating the average 
values of the particle's position and momentum makes 
sense only when this calculation is performed separately 
for the subensembles of transmitted and reflected parti- 
cles. As regards the average values of these quantities cal- 
culated over the whole ensemble of particles, they behave 
non-causally in the course of the scattering process and 
hence cannot be interpreted as the expectation values of 
the position and momentum of a scattering particle. 

By our approach, introducing an observable to describe 
tunneling, without distinguishing transmission and re- 
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flection, is meaningless. Quantum mechanics does not 
imply computing the expectation values of physical quan- 
tities for a tunneling particle, which would be common 
for transmission and reflection. Figuratively speaking, 
it is merely impossible to "pack up" the properties of 
the alternative processes into one common characteris- 
tics. This concerns entirely the dwell time. 

Thus, by our approach, introducing this time scale is 
questionable from the most beginning. Besides, it is ev- 
ident that a proper definition of the transmission time 
for any spatial interval should be valid irrespective of the 
displacement of this interval on the OX-axis. However, 
the definition of the dwell time is evident to violate this 
requirement. Indeed, as is known, this time is defined 
as the ratio of the probability to find a particle within 
the barrier region to the incident flux. We consider this 
definition as a purely speculative one. For the used here 
normalization by the incident flux has no solid physical 
basis. 

Let us assume, for example, that we study the mo- 
tion of a particle in the spatial region of the same width 
but shifted, with respect to the barrier region, toward 
the transmission domain. What flux should be used in 
defining the dwell time for this region? Of course, if 
this spatial region lies entirely in the transmission region, 
then it is naturally to use for this purpose the transmit- 
ted flux (and, as a consequence, a resulting dwell time 
will describe transmitted particles only). However, if the 
shifted spatial region coincides partly with the barrier re- 
gion, then neither incident nor transmitted flux cannot 
be used in the above definition. 

Consider another example. Let the spatial region in- 
vestigated be the right half of the localization region of 
the rectangular potential barrier. At the first glance, the 
dwell time for this case should be the ratio of the proba- 
bility to find a particle within this region to the incident 
flux. However, by our approach, if the particle impinges 
the barrier from the left, then in the case of reflection it 
never enters the right half of a symmetrical potential bar- 
rier. It is evident to be meaningless to take into account 
the reflected part of the quantum ensemble of particles, 
in timing a particle in this spatial region. 



III. RECTANGULAR POTENTIAL BARRIERS 

Now we proceed to a formal comparison of the prop- 
erties of our time scale with those of the phase and dwell 
times. Note flrstly that in any approach the tunneling 
time for the particle with a given value of k can be written 
as mDeff{k)/hk. In this paper we analyze the behaviour 
of Deff{k), in the case of rectangular potential barriers, 
for the phase {Tphase), dwell [Tdweii) and tunneling time 
{Ttun) introduced in our previous paper Q. Of course, 
unlike Tphase and Tf„„, Tdweii describes all particle of a 
quantum ensemble. However, one has to bear in mind 
that in the case of rectangular potential barriers the tun- 
neling times for transmission and reflection are equal. It 



is naturally to expect that in this case the dwell time 
should coincide by value with the transmission time. 

Taking into account the expression for Tphase, Tdweii 
(see, for example, and rt„„, one can easily to obtain 
the expressions for the corresponding effective barrier's 
widths. Let 

jn in 

Tphase{k) = jr^Dpiiase{k) , Tdwell{k) — J-^Ddwell{k) , 



and 



Ttunik) = —deffik). 



For the below-barrier case {E < Vq) we have 



D 



phase 



Ddwell — 



1 2ndk'^{n'^ - P) + nf, sinh(2Kd) 

fc2 2Kd{K'^ - P) + kI sinh(2Kd) 
K Ak'^K^ + K.QSmh^{K.d) 



4 [fc^ + kI sinh^ {Kd/2)] [k^ sinh(Kd) - k'^nd] 
K Ak^ k'^ + K'^sinh^ (Kd) 



where k — ^2m(Vb — E)/h^. 

For the above-barrier case [E >Vo) 



_ 1 2KdP{K^ + fc2) - sin(2Kd) 
K 4fc2K2 + k4 sin^(K(i) 

k^ 2Kd{K^ + k^) - kI sin(2Kd) 



Dd 



K Ak^ k'^ + nfj sin^ (nd) 



^eff 



4 [k"^ - (inl sin^ {nd/2)\ [Pud - finl sin(Kd)] 



4fc2K2 -I- ^4 sin2(K(;) 



where k y/2m{E — Vq)/^?; /3 = 1 if Vb > 0, otherwise, 
/3 = — 1. In both the cases kq = \/2m\VQ\/h?. 

In 01 we treated the case when the potential barrier is 
localized in the region [a, 5], and the CM of the incident 
wave packet to describe all particles is, at i = 0, at the 
point a; = 0; a > 0. Remind that in accordance with the 
standard wave-packet analysis namely this spatial point 
should be considered as the average starting point both 
for transmitted and for reflected particles. However, by 
our approach, this is not the case. As was shown in '/|, 
the average starting points, x%^^i and x[.ll^i, of transmit- 
ted and reflected particles, respectively, do not coincide 
with that for all particles. 

For the above initial condition, in the case of symmet- 
rical potential barriers, we have a;*jQ^t(fc) = x'Jllj.^{k) — 
Xstart{k) where 



Xstart{k) — 2 



Kq (k^ — k^) sinh(Kd) -|- fc^Kdcosh(Kd) 
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Xstart{k) = -2/3 



Kg (k + k ) sin(K(i) — k Kdcos{Kd) 
K Ak'^K'^ + K^sin'^{Kd) 



ior E < Vq and E > Vq, respectively. Note, these quan- 
tities are such that (see 0) 



(1) 



from which a formal connection between the phase-time 
concept and ours is seen explicitly. 



IV. CHARACTERISTIC TIMES IN THE LIMIT 
OF WEAK SCATTERING POTENTIALS 

As is known, by the standard in quantum mechan- 
ics timing procedure (based on timing the CM of a 
wave packet), for a free particle with the well defined 
momentum hk the average time spent by the particle 
in the spatial region of width d is equal to Tfree{k\d) 
where Tfree{k',d) = md/hk; m is the particle's mass, 
k = \p2jmEjW ^ E is the particle's energy. This expres- 
sion is valid for any value of k. In particular, in the limit 

— > 0, Tfree{k; d) diverges as fc~^. 

It is obvious that when the potential energy of a 
particle in the barrier region diminishes (this case is 
named here as the limit of weak scattering potentials) 
then Deff{k) should approach d, for any value of fc. In 
the limit of weak scattering potentials the fc-dependence 
of a true tunneling time must approach the function 
Tfreeik; d). In particular, for infinitesimal potentials the 
transmission time should diverge as k~^ when k —> 0. 
This requirement should be considered as a touchstone 
in solving the tunneling time problem. For example, in 
the case of the rectangular barrier of width d and height 
Vq this should take place when Vod ^ 0. 

At the first glance, all three times behave properly in 
the limit of weak scattering potentials. Indeed, in all the 
cases, if fc 7^ 0, Deff{k)/d 1 when Kgd 0. But, as it 
was stressed above, this should be valid also for however 
small values of k. As will be seen from the following, only 
Tfmi obeys this requirement. 

One can easily show that for Atpd ^ and fc = 



D. 



phase 



D. 



dwell 



Koc'tanh(Koc^) 



0, 



eff 



■ tanh 



d Kod \ 2 y ' d Kodsmh{nQd) 
for Vo > 0; for Vq < 

phase 2 DdyjQii 



(2) 



(3) 



Ko(itan(Koc?) ' d 



0, (4) 



deff _ _2_ / Kod\ Xstart _ 2 

d Kod \ 2 y ' d KQdam(Kod) 



(5) 



Note, in the long- wave limit, D dwell Id = irrespective 
of Kgd. The corresponding limits for the phase time and 
Ttun depend on n^d. Moreover, for fc = 0, as K^d — s- 0, 
D phase /d CO for barriers (Vb > 0), and D phase /d 
— oo for wells (Vq < 0). As regards our definition, in this 
limit deff /d = I both for barriers and wells. 

Figs. 1-4 show the dependence of Deff/d on E/Vq, 
for all three characteristic times, for weak scattering po- 
tentials. Figs. 1, 2 correspond to the narrow {d = 
0.5nm) rectangular barrier (Vb = 0.25eV) and well (Vb = 
~Q.25eV), respectively. Figs. 3, 4 correspond to the 
wide {d — 5Qnm) rectangular barrier (Vb = 0.00025ey) 
and well (Vb = — 0.00025el/), respectively. In all cases 
TO = 0.067TOe where me is the mass of an electron. 

So, in the limit of weak scattering potentials deff/d = 
1 for all values of k. That is, our definition of the tunnel- 
ing time guarantees the passage to the free particle case, 
in this limit. However, this is not the case for the phase 
and dwell times. By these concepts, in the limit of weak 
scattering potentials the time spent by a slow particle in 
the barrier region differ essentially from Tfree{k', d). In 
particular, in comparison with Tfree the dwell-time con- 
cept predicts anomalously short times spent by a slow 
particle in the barrier region. 

To explain this fact, let us remember once more that 
the dwell time was introduced as the ratio of the probabil- 
ity to find a particle in the barrier region to the incident 
probability flux. It is clear that for kq ^ 0, in the long- 
wave limit, the number of particles in the barrier region 
is proportional to k^. That is, a particle with a however 
small value of k does not enter the barrier. Taking also 
into account the fact that the incident flux '--^ A: in this 
limit, we obtain Tdweii ~ k instead of Tdweii ^ k~^. 

This property of the dwell time is kept for a however 
small value of Kod. That is, strictly speaking, the dwell- 
time concept does not imply the passage to the case of 
a free particle, in the limit of weak scattering potentials. 
This fact evidences that the dwell time is ill-defined and 
cannot serve as the characteristic time to describe the 
dynamical properties of a tunneling particle. 

As is seen from the figures, the phase-time concept does 
not guarantee the above passage, too. For slow particles, 
\Dphase\/d may be however large, being negative by value 
in the case of wells. 

Note, to compare the properties of the phase times and 
ours is of a particular importance. "Where do transmit- 
ted particles start from, on the average?" is the main 
intriguing question to arise in this case. Remind, in con- 
trast with the phase-time concept to imply that in the 
above setting the tunneling problem transmitted parti- 
cles start, on the average, from the point x = 0, our 
formalism says that this point is x start- 

Exps. (O and show that in the limit of weak scat- 
tering potentials the phase-time concept predicts an ab- 
normal divergence of the transmission time at A: — *■ 0. 
We have to stress that such a behaviour of the phase 
time takes place even if the transmission coefficient, T, 
approaches unit. Indeed, let us consider the case when 
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kq ^ and k = d/\^ ^ 0; A is fixed. One can easily 
show that in this case 



T = 



that is, T is constant in this hmit. For Vo > we have 



D 



dwell 



D 



phase 



1 

d2' 



for Vo < 0, 
Dd 



well 



= 0, 



D 



phase 



1 



In both the cases 



1, 



start 



-D 



phase- 



As is seen, if X^n^ ^ 4, the barrier is transparent in this 
limit. What is more important, for wells Ddweii = even 
when T « 1. As regards D phase, it diverges in this limit, 
both for barriers and wells. This means, in turn, that in 
this case Xgtart in our approach diverges too (of course, 
for finite wave packets the average value of Ix^tartl does 
not exceed the wave-packet's half- width). 

Figs. 5-8 show the time evolution of the wave pack- 
ets full, ^tr and ^ref to describe all, transmitted 
and reflected particles, respectively. The potential bar- 
rier considered is a rectangular well: Vq — ~712eV, 
d = 1.08 X 10~^nm; in all cases m = 0.067me. At t — Q 
(see Fig. 5) the state of a particle is described by the 
Gaussian wave packet whose half-width equals 15nrn, the 
average particle's energy equals 0.00641eT^; a = 70nm. 
In this case the norm, R, of equals 6.5 • 10~^. 

As is seen from Figs. 6-8, this weak scattering potential 
unexpectedly strongly changes the shape of the "full" 
wave packet. When the peak to correspond, at early 
times, to the CM of the incident wave packet arrives at 
the point lying at some distance of the barrier, a new 
peak to correspond to the CM of the transmitted wave 
packet simultaneously appears behind the barrier region 
(see Fig. 7). 

By the SWPA, the transmission time is negative in this 
case. This property of fuii{x) is usually interpreted as 
the evidence of an ultrafast or even superluminal propa- 
gation of a particle through the barrier. However, by our 
approach, such interpretation of the wave-packet tunnel- 
ing is wrong. The above behaviour of the wave-packet's 
peaks is inherent to interference maxima. It is evident 
that the expectation value of the particle's position can- 
not behave in such a manner. In our approach, for any t, 
this value for transmitted particles should be calculated 
over the wave packet "iltr whose motion near the barrier 
region is regular. 

The above property of the "full" wave packet shows 
once more that there is no direct causal link between the 
transmitted and incident wave packets. By our approach. 



the tunneling process is a combined stochastic process 
consisting from two alternative ones, transmission and 
reflection, provided that fuii — ^tr + ^re/, where 
is the wave function to describe transmission and is 
that to describe reflection (see 0)- 

It is important to remind here that the norms of 5* full 
and '^tr are practically equal in the case investigated. 
However, due to the interference terms (which are pro- 
portional to Vr, rather than R) the difference between 
the behaviour of the elementary state "^tr and combined 
state full (see '7]) is essential. In particular, unlike the 
"full" wave packet, the shape of |^'ir(a;)P remains prac- 
tically unaltered during the scattering event (see Figs. 
6-8). The distance (« 2nm) between the average start- 
ing points of all and transmitted particles is well larger 
than the width of this transparent barrier. 

Remind (see ) , the incoming wave to describe trans- 
mitting particles with a given value of k differs from that 
for all particles by the factor exp(±i7(fc)) where 7(fc) = 

= |i?7[2v^]| 



arctan -yAR/T; R = 1 — T. Since |xstart| 
this shift of the average starting point may be very large 
even if the reflection coefficient R{k) is small by value; for 
the key role here is played by R'{k); the prime denotes 
the derivative on A;. A similar situation arises in the case 
of a resonant tunneling. 



V. RESONANT TUNNELING 

One can easily show that the dwell and phase times, 
distinguished drastically in the long- wave limit, coin- 
cide with each other in the case of a resonant tunneling. 
Namely, for E > Vq, for Kd = mr [n — 1,2,.. .), where 
k = \J /?Kq + v?n'^ /d^, we have 

Dphase D dwell -\- k'^ ^ ^ (jK^Cp 



d d 2k2 

At the same time, if n is even, then 



= 1 



2n27r2 ■ 



d n^TT^ ' 



if n is odd, 



besides, 



start 



(-ly 



2n27r2 



Note, near the resonance point kr the transmission coeffi- 
cient can be written as T{k) = [1 -|-ao(A; — fc^)^]"^, where 
ao is a length to characterize the resonance. It is evident 
that IxstartI =0,0- So that, the narrower the resonance 
peak on T{k), the larger is the value of l^startl- 
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VI. ABOUT THE POSSIBILITY OF 
MEASURING THE AVERAGE STARTING 
POINT OF TRANSMITTED PARTICLES 

Of course, a new and important result to arise in the 
framework of the separate description of transmission 
and reflection is that the average starting point of trans- 
mitted (and reflected) particles does not coincide with 
the initial position of the CM of the incident wave packet 
to describe the whole ensemble of particles. From the 
theoretical point of view, we deal here with three wave 
fields, "iftr and ^ref, where "i/ fuii = ^'tr + *re/- 

Thus, it is not surprising that due to the interference be- 
tween \l/tr and V^re/ the starting points of the CMs of 
the last two wave packets do not coincide, in the general 
case, with that of the incident wave packet. 

This means, in particular, that the average starting 
point of transmitted (or reflected) particles and that of all 
particles cannot be measured simultaneously. Let Oinc 
be an observer to study particles before the scattering 
event. Besides, let Otr and Oref be observers to study 
transmitted and reflected particles, respectively. It is ev- 
ident that from the viewpoint of Omc particles start, on 
the average, from the point a; = 0. This observer cannot 
separate transmission and reflection. On the contrary, 
Otr (Oref) cannot measure the average starting point 
for all particles, but he can measure the average starting 
point of transmitted (reflected) particles. As is follows 
from our approach, he must find that particles transmit- 
ted by the barrier start, on the average, from the point 

start • 

Of course, in verifying our formalism, measuring the 
value of Xstart is of great importance. As in we will 
exploit for this purpose the Larmor precession of a spin- 
1/2 particle in a small magnetic field. However, our aim is 
directly opposite, because timing the scattering particle 
have already been performed in our approach. 

Let an infinitesimal constant uniform magnetic field B 
be applied along the z-axis, everywhere outside the inter- 
val [a — l,b + I] on the x-axis; we assume that a — I > 0, 
moreover a — I ^ Iq and I ^ Iq; Iq is the half- width 
of the incident wave packet. When a spin- 1/2 parti- 
cle moves in this region, the axis of its (average) spin 
will rotate around the z-axis with a constant frequency 
lul ~ g^B/h, where g is the gyromagnetic ratio, fj, is 
the absolute value of the magnetic moment of the par- 



ticle. Of course, we assume that at t = the spin axis 
is not parallel to the z-axis. Moreover, we assume also 
that at this moment the expectation values < Sx > and 
< iSj, > of the X- and y-components of the spin are equal: 



< Si"^ > and < S^"' >. 

Let the spin of a particle with the well-defined value of 
k (the corresponding wave packet is narrow in fc-space) be 
detected at that instant of time, tdet, when the CM of the 
transmitted wave packet arrives at the point x — b + L; 
L ~ I ^ Iq. Thus, on the average, the particle moves 
under the magnetic field during the time Ati„ -I- Atout, 
where AUn = ^{a - I - xstart), Atout = ^{L - I). 
Hence, the expectation values < > and < Sy > for 
transmitted particles, at the instant of time t = tdet, read 
as 



(0) 



V2 



■ cos 



mojL 



hk 



{a + L-2l 



< Sy ^ > 

< Sy >= ^ — sm 



V2 



——(a + L~2l- Xstart) + T 

hk 4 



From this it follows that for this moment in the case of 
equal values of < S^^ > and < Sy^^ > we have 



Xstart — a + L — 21 



hk 



TT f< Sy> 

arctan — 

4 V < > 



(6) 



Note, according to the SWPA the expression in the right- 
hand side of © should be equal to zero. 

Thus, Exp. © can serve as the basis for the experi- 
mental checking of our approach. They provide the way 
of calculating the value of Xstart from experimental data 
on measuring the particle's spin. We have to stress once 
more that these expressions describe the case when an 
external magnetic field is applied to the spatial regions 
where the incident and transmitted wave packets evolve 
freely. 

Note, Exp. ® is valid also for finite wave packets. 
However, now the expectation values of k for the inci- 
dent (< k >inc), transmitted (< k >tr) and reflected 
(< k >ref) wave packets are different. Thus, fc in © 
should be replaced in this case with < k >tr- 
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Figure captions 

Fig. 1. Df,ff/d versus E/Vq for the narrow [d 
Q.bnm) rectangular barrier (Vq = 0.25ey). 



Fig. 2. Dgff/d versus E/Vq for the narrow [d = 0.5nm) 
rectangular well (Vq = —0.25eV). 



Fig. 5. The x-dependencc of (solid line) which 

represents the Gaussian wave packet with Iq = 15nm and 
the average kinetic particle's energy 0.00641ey, as well 
as |*trP (open circles) and |^'re/P (dashed line) for the 



rectangular well (Vb 
a = lOnm); t = 0. 



-7UeV, d = 1.08 x IQ-^nm, 



Fig. 6. The same as in Fig. 5, but t = 29ps. 



Fig. 3. Deff/d versus E/Vq for the wide {d = 50nm) 
rectangular barrier (Vb = 0.00025eV). 



Fig. 7. The same as in Fig. 5, but t = 33.5ps. 



Fig. 4. Deff/d versus E/Vq for the wide {d = 50nm) 
rectangular weU (Vb = -0.00025ey). 



Fig. 8. The same as in Fig. 5, but t = 38ps. 
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